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Abstract 

The physics of the initial conditions of heavy ion collisions is dominated by the nonlinear gluonic 
interactions of QCD. These lead to the concepts of parton saturation and the Color Glass Conden- 
sate (CGC). We discuss recent progress in calculating multi-gluon correlations in this framework, 
prompted by the observation that these correlations are in fact easier to compute in a dense system 
(nucleus-nucleus) than a dilute one (proton-proton). 
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1. Introduction 

Bulk particle production in relativistic collisions around midrapidity originates from small x 
degrees of freedom, predominantly gluons, in the wavefunctions of the colliding hadrons or nuclei. 
At large energies these gluons form a dense system characterized by a saturation scale Q s . The 
degrees of freedom with px < Q s are fully nonlinear Yang-Mills fields with large field strength 

~ 1/g and occupation numbers ~ l/a s ; they can therefore be understood as classical fields 
radiated from the large x partons. Because of their large longitudinal momentum, the large x 
degrees of freedom are effectively "frozen" during the interaction. They can be described as random 
color charges drawn from a classical probability distribution Wy[p] that depends on the rapidity 
cutoff y = In 1/x separating the large and small x degrees of freedom. The dependence of W y [p] on 
y is described by a Wilsonian renormalization group equation known by the acronym JIMWLK. 
Note that while this description is inherently nonperturbative, it is still based on a weak coupling 
argument, because the classical approximation requires a s (Qs) to be small and therefore Q s 3> 
A QCD . The Color Glass Condensate (CGC, for reviews see [l|]) is a systematic effective theory 
(effective because the large x part of the wavefunction is integrated out) description of the classical 
small x degrees of freedom. 

The term glasma 0] refers to the coherent, classical field configuration resulting from the 
collision of two such objects CGC. The glasma fields are initially longitudinal, whence the "glasma 
flux tube" 0, 0] picture. More importantly for computing multigluon correlations, they are boost 
invariant (to leading order in the QCD coupling) and depend on the transverse coordinate with 
a characteristic correlation length 1/Q S . There are several signals in the RHIC data 0, 0] that 
point to strong correlations originating from the initial stage of the collision. The glasma fields 
provide a natural framework for understanding these effects, although much work is still left to do 
in understanding the interplay with purely geometrical effects from the fluctuating positions of the 
nucleons in the colliding nuclei 0, 0] • 

We shall first describe some general observations on computing multigluon correlations in the 
glasma, arguing in Sec. [2] that they are in some sense simpler to compute in a collision of two 
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Figure 1: Left: Building block, Lipatov vertex coupled to two sources. Right: combinatorics of the sources. The 
combinatorial problem is to connect the dots on the upper and lower side (left- and right moving sources) pairwise. 



dense, saturated nuclear wavefunctions than in the dilute limit (see Ref. [9( for a more formal 
discussion). We shall then, in Sec. [3] discuss one application of these ideas to computing the 
multiplicity distribution of gluons in the collision before moving to the leading In 1/x rapidity 
dependence of the correlation in Sec. [H 



2. Multigluon correlations in the glasma 

The gluon fields in the glasma are nonperturbatively strong, ~ 1/g. This means that the 
gluon multiplicity is N ~ l/a s - For a fixed configuration of the classical color sources it is well 
known that the multiplicity distribution of produced gluons is Poissonian, i.e. (iV 2 ) — (iV) 2 = (iV). 
In this case the correlations and fluctuations in the gluon multiplicity are all quantum effects that 
appear only starting from the one-loop level, i.e. suppressed by a power of the coupling constant a s . 
The computation in the CGC framework does not end here, however. To calculate the moments of 
the gluon multiplicity distribution one must first calculate the gluon spectra for fixed configuration 
of the color charges p and then average over the probability distribution W^[p(a;_|_)]. For the nth 
moment of the multiplicity distribution, i.e. an n-gluon correlation, the leading order result is 
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where the subsrcipt "LO" refers to the single gluon spectrum evaluated from the classical field 
configuration corresponding to a fixed configuration of color charges. This averaging, even after 
the subsequent subtraction of the appropriate disconnected contributions, introduces a correlation 
already at the leading order in a s , i.e. enhanced by an additional l/a s compared to the quantum 
correlations. A natural example is the negative binomial distribution that we shall discuss below, 
whose variance is (N 2 ) — (-/V) 2 = (N) 2 /k + (N). One must emphasize here that although these 
contributions arise as formally classical correlations in the effective theory that is the CGC, they 
are physically also quantum effects, where the weak coupling is compensated by a large logarithm 
of the energy that has been resummed into the probability distribution W y [p(x±)]. In this sense 
the leading correlations are present already in the wavefunctions of the colliding objects. 



3. Multiplicity distribution 

We can then apply this formalism to the calculation of the probability distribution of the 
number of gluons in the glasma 10|. We shall assume the "AA" power counting of sources that 
are parametrically strong in g, but nevertheless work to the lowest nontrivial order in the color 
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Figure 2: Numerical evaluation of the two-gluon correlation strength in the MV model as a function of the two 
transverse momenta p± and q±. Left: the correlation strength on the "near side" \tp p — tp q \ < n/2, right: the "away 
side" \tp p — <p q \ > 7r/2. 



sources. Formally this would correspond to a power counting p ~ g 6 ^ 1 with a small e > 0. In this 
limit, as we have discussed, the dominant contributions to multiparticle correlations come from 
diagrams that are disconnected for fixed sources and become connected only after averaging over 
the color charge configurations. The corresponding two gluon correlation function was computed 
in Ref. 0] and generalized to a three gluons in Ref. 11]. We shall here sketch the derivation [lF 
of the general n-gluon correlation in this simplified limit. 

Working with the MV model Gaussian probability distribution 



W[p] = exp 



- p a (x ± )p a (x ± ) 
d. x± t-= 



(2) 



computing the correlations and the multiplicity distribution in the linearized approximation is a 
simple combinatorial problem. Each gluon is produced from two Lipatov vertices (see fig. [I] left), 
one in the amplitude and the other in the complex conjugate. The combinatorial factor is obtained 
by counting the different ways of contracting the sources pairwise (see fig. [T] right). The dominant 
contributions are the ones that have, in the dilute limit, the strongest infrared divergence which is 
regulated by the transverse correlation scale of the problem, Q s . When integrated over the momenta 
of the produced gluons one obtains the factorial moments of the multiplicity, which define the whole 
probability distribution. It can be expressed in terms of two parameters, the mean multiplicity 
n, and a parameter k describing the width of the distribution. The result of the combinatorial 
exercise is that the qih factorial moment m q (defined as (N q ) minus the corresponding disconnected 
contributions) is 
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Here S± is the transverse area of the system These moments define a negative binomial distribution 
with parameters k and n, known as a phenomenological observation in high energy hadron and 
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Figure 3: Left: Comparison of a two-particle correlation computed using eq. ^ supplemented with a short- 
range correlation contribution from PYTHIA with PHOBOS data. Right: Rapidity correlation at LHC energies 
fcx-factorization approximation. Plots fromfljl]. 



nuclear collisions already for a long time. Also the numerical magnitude of the parameter k obtained 
from the saturation scale agrees very well with both pp and AA data [r3 |. 

This calculation predicts that the k parameter should increase with energy, unlike what is seen 
in UA5 data. Multiplicity fluctuations in proton-proton collisions at lower energies are still mostly 
dominated by the dilute edge of the collision system, causing a Poissonian nature (i.e. k — > oo) 
of low energy particle emission. Our calculation formally assumes Q^Sj_ 1, and we expect the 
growing behaior of k with energy to eventually take over at high enough energy. Some signs of this 



are already visible in the LHC yfs = 7 TeV data[13l]. 



In terms of the glasma flux tube picture this result has a natural interpretation. The transverse 
area of a typical flux tube is l/Qg, and thus there are Q^S± = A^ft independent ones. Each 
of these radiates particles independently into N c 2 — 1 color states in a Bose-Einstein distribution 
(see e.g. [3). A sum of k w Nft(N c 2 — 1) independent Bose-Einstein-distributions is precisely 
equivalent to a negative binomial distribution with parameter k. A numerical evaluation [l5| of 
the second moment of the distribution, parametrized in terms of 

t ^ n2c ( d * N dN dN \ I dN dN ^^ 

k2(p±,q±) = Q s s± (^dV " d^dV J / d^IdV (4) 



is shown in Fig. [2j It confirms the expectations of [10| that this ratio is of order one and depends 
only weakly on the momenta , . 



4. Rapidity dependence 

The general discussion of Sec. [2] on the different nature of multigluon correlations in the "A A" 
case applies also to the rapidity dependence. Until now we have only been discussing gluon pro- 
duction in a rapidity interval smaller than l/a s . For this we needed only the correlations between 
the color charges p(x±) measured at this same rapidity. To understand the rapidity dependence 
of the correlations one needs also the correlation between color charges at different rapidities, 
{Py( x ±)Py'(y±)) ■ Also this information is contained in the JIMWLK renormalization group evolu- 
tion, at least to leading lnl/x accuracy [l7t ] . An intuitive description of the resulting correlations 
is provided by the formulation of JIMWLK as a Langevin equation in the space of Wilson lines 
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formed from the color charges. In this picture the evolution proceeds in individual trajectories 
along an increasing rapidity. A first attempt of a realistic estimate of the rapidity dependence 
of two-gluon correlations is performed in Ref. [3]. Evaluating the two gluon correlation in a di- 
lute limit in a fej_-f actor ized approximation, but keeping the general structure resulting from the 
JIMWLK evolution leads to the following expression: 
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(5) 



Note the very different structure of this correlation compared to one where the gluons would 
be produced from the same diagram for fixed sources. The two gluon correlation function is 
proportional to the product of four unintegrated gluon distributions, with three of them evaluated 
at the rapidity of one of the produced gluons and only one at the other. This structure is a 
direct consequence of the nature of JIMWLK evolution. The resulting correlation is compared to 
PHOBOS data in fig. [3l The /cj_-factorized approximation gives a very inaccurate description of the 
gluon spectrum in the transverse momentum regime p± ~ Q s where the bulk of the particles are 
produced [3]. Equation ([5]) has also been derived in the approximation, true only in the linearized 
case, that the unequal rapidity correlation of two color charge densities is equal to the unintegrated 
gluon distribution at the smaller one of these rapidities. As of yet there is no calculation of how 
much this approximation is violated in the full JIMWLK evolution. The results presented in 
fig. [3] are therefore not the final word on the subject, alhough it is reassuring that such a simple 
approximation seems to agree rather well with the experimental result. 



5. Forward-backward multiplicity correlation 

Let us conclude by a few remarks on the somewhat puzzling STAR data rcl on backward- 
forward charged hadron multiplicity correlations, following the discussion in Ref. [8(. The reported 
forward-backward correlation shows a rapid increase as a function of centrality, and appears to 
have a strength which cannot be explained by a superposition of pp interactions. The measured 
quantity is usually reported as a correlation coefficient between the charged hadron multiplicities 
measured in a forward and a backward rapidity bin (N F and N B ) within the STAR TPC coverage 
of two units in pseudorapidity: 

(N F N B ) - (N F )(N B ) 

(N F 2 )-(N F ) 2 ' W 

where due to the symmetrical placement of the bins around midrapidity the variance {N F 2 ) — (N F ) 2 
is the same for both rapidity bins. If the expectation value here were taken over all the events in 
a centrality bin, correlations would be generated by the different impact parameters (or numbers 
of wounded nucleons) possible within such a This is, however, a too simplified picture of 



One expects the charged multiplicity to be strongly correlated with impact parameter and if the impact parameter 
itself can have significant variation within a fixed centrality bin, then spurious correlations whose only origin is the 
geometry of the collision would be generated. 
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Figure 4: Simple parametrization of a measured forward-backward correlation at fixed reference multiplicity. Left: 
dependence on strength of the long range correlation K is minimal. Right in stead, the observed correlation strength 
mostly depends on the short range correlation a, i.e. the fluctuations of the multiplicities that are uncorrelated 
between different rapidities. 



what is actually done in the analysis. The nature of the actual quantity reported by STAR is not 
very explicitly described in the paper; to understand it one must decipher the inocuous-looking 
statement @] ". . . (N F , N B , N F 2 , N F N B ) was obtained on an event-by-event basis as a function of the 
event multiplicity [JVr] . . . " . In other words the correlation coefficient is not measured averaging 
over all the events in the bin, but over events with a fixed multiplicity in a third, reference, 
rapidity window, i.e. 

b = (N F N B ) NK -(N F ) Nli (N B ) NK 

(N F 2 ) Nn - (N F ) 2 Nn 

Thus the measured quantity does not describe the relation between two but three correlated mul- 
tiplicities. For the largest separations between the forward and backward rapidity windows the 
reference window is between the two, so there is no reason to assume it to be less correlated with 
the F, B windows than these are with each other. The analysis in Ref. [8y shows that in the limit 
of a maximal correlation between multiplicities at different rapidities the coefficient "6" of eq. ([7J) 
reaches a maximum value of 1/2. The results of a more detailed parametrization including the 
fluctuations in centrality and two parameters describing the magnitudes of the uncorrelated short 
range fluctuations ~ a and the long range correlation ~ K are shown in fig. [5j The result remains 
that while the measured correlation never exceeds 1/2, its magnitude depends more on the uncor- 
related fluctuations than the actual long range correlation. The experimental measurements 0] 
saturate and, for the most central bins, even exceed this limit derived on very general grounds. 
Although they do seem to point to very strong correlations, this inconsistency makes their in- 
terpretation in terms of any microscopic origin, the glasma or something else, of the correlations 
difficult. 



6. Conclusion 

Most experimental observables do not probe the glasma initial state of directly, because the 
system goes through a complicated time evolution before the hadronization stage. A good candidate 
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for an experimental probe giving direct access to the initial state is provided by different kinds of 
correlation measurements. These have indeed been a focus of both experimental and theoretical 
activity recently. We have argued here that the glasma picture of the initial stages of a heavy ion 
collision is the natural framework to understand the origin of these correlations. 
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